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(Dated: 2 November 2016) 

Renewables will soon dominate energy production in our electric power system. And 
yet, how to integrate renewable energy into the grid and the market is still a subject 
of major debate. Decentral Smart Grid Control (DSGC) was recently proposed as a 
robust and decentralized approach to balance supply and demand and to guarantee a 
grid operation that is both economically and dynamically feasible. Here, we analyze 
the impact of network topology by assessing the stability of essential network motifs 
using both linear stability analysis and basin volume for delay systems. Our results 
indicate that if frequency measurements are averaged over sufficiently large time in¬ 
tervals, DSGC enhances the stability of extended power grid systems. We further 
investigate whether DSGC supports centralized and/or decentralized power produc¬ 
tion and find it to be applicable to both. However, our results on cycle-like systems 
suggest that DSGC favors systems with decentralized production. Here, lower line 
capacities and lower averaging times are required compared to those with centralized 
production. 

PACS numbers: 05.45.Xt: Oscillators, coupled, 89.75.-k: Complex systems, 84.70.+p: 
High-current and high-voltage technology: power systems; power transmission lines 
and cables, 88.05.Lg: Cost, trends in renewable energy 
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I. INTRODUCTION 


The ongoing climate change is forcing us to shift our power generation from fossil power 
plants towards renewable generation 1 . In the last years, renewable energy technology devel¬ 
opment and policy support led to a tremendous increase in the share of Renewable Energy 
Sources (RES). In 2014, Germany covered 27.8% of its gross electricity consumption with 
RES 2 . Still, large conventional power plants dominate the power grids: transmission lines 
connect large plants with regional consumers in a locally star-like topology. With more 
renewable power sources entering, the grid topologies become more decentralized and more 
recurrent due to the distributed generatiorP. In such a scenario, consumers may act as pro¬ 
ducers and consumers at the same time, so-called prosumer P and electricity transport is no 
longer unidirectional. 

A known challenge of many renewable sources is their volatile naturP 7 . Fluctuations 
occur on different time scales, including seasonal, inter-da^ 8 and intra-second fluctuation^ 9 . 
This requires radical changes in the control and design strategies of electric power grids as 
well as market innovations to ensure cost effectiveness. Therefore, a need for more flexibility 
options for power market supply and balancing energy 10 arises because the fluctuating RES 
cannot guarantee power supply with the certainty conventional plants could. In this regard, 
it is most important to identify options that are both cost efficient and system stabilizing. 
So far, the framework of power market design and power grid stability with its long planning 
horizons does not satisfy the need for sufficient flexibility optionPR 

Different smart grid approaches have been proposed to present ways to match supply 
and demand in such a fluctuating power grid. However, economic and political feasibility 
and market integration are often missed out. A key idea of various smart grid concepts 
is to regulate the consumers’ demand 1 -^, a massive paradigm shift compared to the current 
power grid operation schemes 13 14 . Many proposals for smart grids are based on sufficient 
information and communication technology infrastructure, see, e.g.J® oi 16 . However, such 
a centralized system would raise questions of cyber security and privacy protection 17 18 and 
several studies highlight the cost burden these proposals implicate. 

In contrast, an alternative approach without massive communication between consumers 
and producers directly utilizes the grid frequency to adjust production and consump¬ 
tion. The frequency increases in times of power excess while it decreases in times of 
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underproduction 20 21 . A novel smart grid concept, Decentral Smart Grid Control (DSGC), 
was introduced irP, based on earlier ideas by®l, and its mathematical model proposed and 
analyzed in 23 . Using DSGC prosumers control their momentary demand on the basis of 
the grid frequency which can easily be measured everywhere with cheap equipment. Ref. 23 
demonstrates that DSGC enhances the stability of the power grid dynamics if the frequency 
measurements are averaged over sufficiently long time intervals. Yet, so far, only very small 
networks were investigated. Hence, the impact of grid topologies on power grid dynamics 
with DSGC constitutes a widely open research question. 

Here we analyze the impact of network topology by assessing the stability of essential 
network motifs using both linear stability analysis of delay systems and determining basin 
volume. Furthermore, we address the question, how grid stability changes when generation 
is decentralized. The article is structured as follows. First, we present a dynamical model for 
power grid dynamics and present the concept of Decentral Smart Grid Control (DSGC) 23 
to control a power grid in section |TT[ In section HI we briefly summarize linear stability 
and basin volume measures for such delayed systems. The stability results of DSGC are 


then presented for a star motif in section IV A where we discuss the destabilizing resonance 
and rebound effects and how stable grid operation remains possible. Using linear stability 
analysis, we investigate the effect of decentralized power generation in cyclic and square 


lattice grid motifs in section |IV B| The results suggest that DSGC works successfully for 
centralized as well as decentralized production, where grids with decentralized production 
require lower line capacities than centralized ones. 


II. COUPLED OSCILLATOR MODEL WITH DECENTRAL SMART GRID 
CONTROL 

To model the frequency dynamics of a large-scale power grid, we consider an oscillator 
model based on the physics of coupled synchronous generators and synchronous motors, 
sed 24 30 for details. This model is similar to the ’’classical model’ 3 - and the ’’structure pre¬ 
serving model’ 32 from power engineering. 

The state of each machine i G {1,..., A} is characterized by the rotor angle &i(t) relative 
to the grid reference rotating at D = 2n x 50 Hz or D = 2tt x 60 Hz, respectively, and its 
angular frequency deviation Ui = dd.Jdt from the reference. Each machine is driven by a 
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mechanical power P t (t), which is positive for a generator and negative for a consumer. In 
addition, every machine transmits electric power via the adjacent transmission lines which 
have a coupling strength K,j. This coupling strength expresses the maximal possible power 
that may theoretically be transmitted through the power lines. The dynamics of the machine 
i is then given by the equation of motion as 


d^ 

d t 2 


Pi ~ a~ + ^ Kij sin (0j - 6i) Vi G {1,..., IV}, 


( 1 ) 


3 = 1 

where cq is a damping constant. We neglect ohmic loads which should be small compared 
to shunt admittances 33 for the dynamics we consider. We take the moment of inertia to be 
identical for all machines and hence eliminate such moments of inertia in the equation of 
motion for simplicity of presentation. Equation ([TJ) as well as the upcoming equations ([5]) 
and <§ are discussed in more detail irP. 

Decentral Smart Grid Control (DSGC) is based on Demand Response that aims to stabi¬ 
lize the power system by encouraging consumers to lower their consumption in times of high 
load and low production and increase consumption in times of low load but high production. 
Instead of paying a constant price for electric power, consumers are presented with a linear 
price-frequency relation Pi(^) 


d 6i d 9i 

Pi( dd= pa ~ C1 ' dT 


( 2 ) 


to motivate grid-stabilizing behavior. Although consumer reaction might be very complex, 
we assume a linearized power-price relation Pi(pi) 


Pi(Pi ) ~ Pi + c 2 ■ (p* - pa) 


( 3 ) 


by the consumers close to the stable operational state. Plugging (|2} into (|3} and defining 
7 = Ci ■ c 2 leads to a linear response of consumed and produced mechanical power Pi(t) as 
a function of frequency deviation ddj/df: 


Vi € {1,..., IV}, (4) 

where 7 j is proportional to the price elasticity of each node i, i.e., measures how much a 
producer or consumer is willing to adapt their consumption or production, see alsd- 23 ! In 
general, such an adaptation will not be instantaneous but will be delayed by a certain time r 
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FIG. 1. Using linear relations the power becomes a linear function of the frequency 
deviation d6i/dt. (a): We assume a linear price-frequency relation to motivate consumers to sta¬ 
bilize the grid. For example, if the production is larger than consumption, the power grid frequency 
increases. Hence, decreasing prices should motivate additional consumption, (b): Although con¬ 
sumers might react non-linearly towards price-changes (dark blue), we assume a linear relationship 
(light green) close to the operational frequency fl which corresponds to ddi/dt = 0. 

by a measurement and the following reaction. We can now substitute the function P t (t — r) 
from Q for the fixed value Pi in the uncontrolled system ([TJ) and obtain the equation of 
motion 


N 


d 2 #- df) • , iV r\fl. 

= Pi - OLi-^ + ^KijSm^j - 6i) - 7i-^(t - r) Vi G {1, ..., N}, (5) 

3 = 1 

with DSGC including a delayed power adaptation. In 23 it was already shown that such a 
delayed system poses risks to the stability of the power grid for certain delays r. Hence, 
an extension using frequency measurements averaged over time intervals of lengths T were 
introduced to stabilize the power grid regardless of the specific delay. Such averaging yields 


d% _ Mi 

d t 2 1 ai d t 


iV -in 

+ V" Kij sin (9j - 0i )~tp ~ T)dt' 

7=i 1 Jt ~ 


dt 


( 6 ) 


1 = 1 
N 


dQ r _ „ x % 


= Pi - OLi-^ + ^ sin (0j - 9, 

3 = 1 


)-dd{e l {t~r)-e l {t-T- t)) Vi g {i,...,iv}7) 


For what follows, we choose homogeneous averaging time T for all machines, as well as 
similar delays r for all nodes. In addition, we use homogeneous capacities = K for 
all lines to simplify the calculations. In the section IV we apply equations ([5]) and Q to 
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different network topologies and evaluate their stability as a function of the delay r with 
different averaging times T. We hereby treat the averaging time T as a control parameter 
that can be chosen when setting up the system, while the delay r remains as an exogenous 
parameter introduced by the consumers and producers. 


III. ASSESSING ROBUST OPERATION 


Here, we discuss how linear stability analysis and measuring basin volume yield infor¬ 
mation about robust operation, in dependence of delay r and for different averaging times 
T. First, we introduce the fixed point of the system, then discuss linear stability analysis 
of delayed systems and finally point out difficulties when assessing basin volume of a power 
grid system with delay. 

a. Fixed points. To study the stability and the role of the system parameters, we 
analyze the dynamical stability around the steady-state operation of the grid given by the 
fixed point 

= (»;.v;e{i,...,v}, (») 

as obtained by solving 


Since u>i = ddi/dt we directly obtain u ;* = 0 for all i. Hence, we only need to determine the 
angles 6*. A fixed point can only exist, if the grid has a sufficient transmission capacity K t] 
to transmit the power from the producers to the consumer^. The minimal K tJ for which a 
stable fixed point exists is called critical coupling 26 . 

b. Linear stability. Linear stability of a dynamical system is determined by the eigen¬ 
values of its characteristic equation. For systems without delay this is a polynomial obtained 
from the Jacobian of the system but for a delayed system it becomes a quasi-polynomial 
with infinitely many solutions 35 36 . We obtain the characteristic equation by calculating the 
Jacobian of both the non-delayed system, based on equation ([5]) with r = 0, 


J 0 = 


a 

( d a .\ 

a 

( d f) • 

ddi 

\dt U V 

duii 

Id t U X 

a 

(d, t 

a 

( d , , 

ddi 

Id t^o) 

doJi 

Id t U 3 


e M 


2N x 2N 

1 


( 10 ) 
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and the derivatives for the delayed terms involving r, 


Jr = 


a 

ae T , 

a 

ae T . 


( d 0 A 

a 

( d o) 

Idt'W 

du> T ,i 

IdtJj 

(sA) 

a 

dul T ,i 


r) and 

^r,i ~ 



t)2N x 2N 


( 11 ) 


consider exponentially decaying or growing solutions 35 . The stability eigenvalues A are then 
determined by the solutions of the characteristic equation 


p( A) = det(J 0 + e Ar J T — All) = 0. 


( 12 ) 


For the delayed system with averaging, i.e., equation ([6]), we simply calculate the delayed 
Jacobian for the two delays r and T = T + r. Hence, the characteristic equation is given by 


p{ A) = det(J 0 + e Ar J T + e X1 Jf — All ) = 0. 


-AT 


(13) 


We obtain the symbolic expression for the characteristic equation using Mathematical 7 which 
is then also used to numerically determine roots of the characteristic equation, via Newton’s 


method. Equations (12) and (13) have infinitely many solutions but only a finite number of 


those can have a positive real part and those determine the instability of the systemP 8 . 

Our method of finding these eigenvalues works as follows: We start at an arbitrary 
delay r = r samp ii n g > 0 and let Mathematica find approximately 10,000 roots by choosing 
random complex initial conditions for Newton’s algorithm. Afterwards, we delete double 
entries. The obtained eigenvalues are taken as the initial conditions for Newton’s algorithm 
for the next larger delay r = T samp i in g + 0.01s. These eigenvalues then serve as the initial 
conditions for the next delay step etc. Similarly, we obtain eigenvalues for smaller delays 
like r = T samp ii n g — 0.01s by using again the eigenvalues from r samp ii ng as initial conditions. 

Linear stability analysis quantifies whether a fixed point is stable to small perturbations 
and constitutes a fundamental aspect of stability in dynamical system. Assessing the sta¬ 
bility of the system with respect to larger perturbations requires further analysis. 

c. Basin volume. The global stability of a fixed point of a dynamical system can be 
quantified by the volume of its basin of attraction. An estimate for the basin volume Wasin 
is determined numerically using a Monte Carlo method as the ratio of initial conditions 
converging to a stable operation state to the total number of initial conditions, as proposed 
in 39 . Note that delayed systems are infinite-dimensional 35 and do need an initial function 








instead of a single initial condition. We treat this problem by setting the initial function to 
be identical to the initial condition for all times smaller than zero, i.e, 

0(t < 0) = 9(t = 0), oj(t < 0) =u(t = 0). (14) 

Thereby, we can effectively choose initial conditions as they completely define the initial 
function. In the following we take M = 1000 randomly chosen initial conditions into account 
in order to estimate the basin volume’s dependence on the delay time r. 

We are mainly interested in how fluctuations or disturbances in the energy generation 
will influence the system’s dynamics. Hence, we first perturb the producer’s node phase 
angle and angular velocity around its component of the fixed point (see fig. for the 
network topology). In the next series of simulations, we perturb one randomly selected 
customer node around its component of the fixed point. Perturbations are uniformly chosen 
at random from the intervals A 6 1 ,, G [—7r, 7r] and Ac u* G [—30, 30] Hz for the initial angles 
and initial frequencies respectively, similar td 29 40 . We run the simulations for a simulation 
time of f S i m = 1500s. These long and computationally costly simulation times are necessary 
because we observed that for specific values of delay time r, e.g., r = 1.4s, perturbations 
may decay relatively quickly toward the fixed point but later still escalate. 

Note that we only consider so-called single node basin volumes, i.e., we only perturb 
the component the a fixed point of one node. In theory, all nodes could be perturbed 
simultaneously which results in a more complete sample of the phase space. Unfortunately, 
the total phase space volume grows exponentially with the number of nodes, making it 
infeasible to sample the full phase space. 

IV. RESULTS 

We now present results about networks with Decentral Smart Grid Control. First, we 
present and compare the results of linear stability and basin volume analysis of a four node 
star motif. This motif constitutes one of the main building blocks of power grids, since, 
in principle, its effective topology locally resembles a star, the central node being a large 
power plant that supplies the regional consumers in its vicinity 6 26 . Hereby, we discuss 
the destabilizing effects of resonances and the ’’rebound effect” for large delays. Using basin 
volume we present how intermediate delays r benefit the stability. Finally, we consider larger 
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networks and demonstrate how decentralization enhances stability. The parameters of the 
swing equation are calculated from standard literature values 31 10 . In current (European) 
power plants the initial delays have to be smaller than 2 seconds according to European 
regulation^, in practice they will be significantly smaller. However, in future power grids 
additional communication delayd 4 ^ of the order of several hundred milliseconds might arise 
in addition to unknown delays caused by demand response and additional power electronics. 
Hence, we consider a large range of potential delays r E (0, 5)s looking for the boundary of 
acceptable delays. 


A. Stability of the star motif 



a) 
a: 


Eigenvalues of four node network 


o.io 



2 4 6 

delay r [s] 


T=0 s 
T=2 s 


FIG. 2. Resonances and large delays r destabilize the four node system, (a): An 
elementary building block in a power grid with centralized production is shaped like a star. Shown 
is a motif whose linear stability and basin volume we study. The network is formed of one producer 
(green) in the center with power P pro ducer = 3/s 2 and three consumers (red) with power /Consumer = 
— 1/s 2 each, (b): Plotted are the eigenvalues with the largest real part as functions of delay r. For 
no averaging (dark blue curve), stable and unstable regions exist. For an averaging of T = 2s, the 
system is stable for all delays below a critical r c ~ 8s. In ([b]) parameters a = 0.1/s, K = 8/s 2 and 
7 = 0.25/s were applied. 

d. Delays induce destabilizing resonances. Networks with star topology (fig [2^,) exhibit 
stability properties that depend crucially on the delay and the averaging applied (fig. If)- 
Without any averaging (fig. [2 ]d dark blue curve), there are delays r for which the fixed point 
is linearly unstable, i.e., there are eigenvalues with a positive real part Re A > 0. Those 
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eigenvalues exhibit a periodic behavior with respect to the delay r. Operating the power grid 
at a delay r for which we find a positive real part, e.g., r ~ Is, is equivalent to resonantly 
driving the power grid away from the fixed point instead of damping it towards stable 
operation. These destabilizing delays are linked to the eigenfrequency of the oscillators in 
the power grid. If the delay is half the eigenoscillation duration, then it increases amplitudes 
of perturbations instead of damping them. This destabilization only occurs for a < 7 because 
the resonant driving has to be larger than the intrinsic damping of the system, see also 23 . 
Introducing sufficiently large averaging times into the control cures these instabilities (fig. 
§> light green curve); the unstable regions vanish for all delays r < 7s. 

e. Rebound effect for large delays. For delays larger than a critical delay r > r c ~ 8.7s 
the system always gets destabilized, i.e., there is an eigenvalue with Re(A) > 0. This 
rebound effect acts on a longer time scale than the intrinsic oscillations of the power grid 
system and originates from an over-reaction of the attempted damping as we explain below. 
The existence of such a rebound effect is independent of averaging T (fig. [2]) . We determine 
the critical delay without averaging r c to be 


t c 


arccos 



+ 


2irn 


y/l 


2 — ex 2 


n E Z. 


(15) 


This result is obtained by the following considerations. We define the sum of all angles as 
£0 := an d obtain its equation of motion by using eq. ^ as 


d^ 

df 2 


£0(t) 


-afzm-ifmt-T). 


(16) 


The characteristic equation of this equation reads 


p (A) = — ol — ye Ar — A = 0, 


(17) 


where we eliminated a zero eigenvalue A = 0 which arises due to the possibility to shift all 
angles by a constant. For r = 0 the eigenvalue A = —a — 7 is negative as ex > 0 and 7 > 0 
and hence the system is stable with respect to the sum £0. For larger delays r > 0 we set 
A = i ■ £ to obtain the delays for which the stability changes. We get 

— ex — 7 e _, ^ r — if = 0 . (18) 
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Applying complex expansion and separating into real and imaginary parts we obtain 

— a — 7 cos(£t) = 0 (19) 

7Sui(£t) — £ = 0. (20) 


These equations can be solved for r and £ to yield the critical delay as in eq. (15). Note that 


a critical delay r c only exists, if the price adaptation is larger than the intrinsic damping 
of the system 7 > a. Following straight-forward calculations we can prove that eigenvalues 


obtained from eq. (17) always destabilize the system, i.e., their real parts are positive for 


all delays larger than the critical one, 


Re(A(r)) >0 Vr > r c 


( 21 ) 


These results hold for all network topologies, since we needed no assumptions regarding the 
coupling matrix Kjj or the power production P l . Predicting the precise scaling of the critical 
delay as a function of the averaging time T is not easily possible but an approximation for 
small £T is obtained as 

a/T 2 7 2 — 4 arctan (a+7 ) A / T 1 — 4 
r c (T) «- 


{2+T~))yJ a 2 —"/ 2 


( 22 ) 


y/a 2 — 7 2 

which is a decreasing function in T for parameters a, 7 ,T > 0. Hence, increasing averaging 
time T causes the rebound effect to occur for smaller delays r. 

We conclude that the delay r has to be smaller than a critical value r c to ensure stability. 
This critical value depends only on the intrinsic damping a and the price adaptation 7 and 
decreases for increasing averaging T, while it is valid for all network topologies. Hence, to 
avoid problems with large delays, we have to enforce all actors of the power grid to react 
within less then this critical delay r c or need to ensure that intrinsic damping is larger 
than the price adaptation: a > 7 . For the next section, we restrict ourself to the interval 
r G [0, 5]s to avoid this destabilizing rebound effect. 

/. Intermediate delays benefit stability. With the help of linear stability we observed 
that delays induce destabilizing resonances which can be suppressed by prosumers respond¬ 
ing to averaged frequency data. At the same time large delays destabilize the system by 
introducing a rebound effect. These results are supplemented by information from basin 
volume analysis. For DSGC with averaging (fig. and c), we demonstrate how interme¬ 
diate delays r are beneficial for the stability of the system. The basin volume increases 
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'/Basin consumer 
'/Basin producer 
Re [A] 


'/Basin consumer 
'/Basin producer 
Re [A] 


'/Basin consumer 
'/Basin producer 
Re [A] 


FIG. 3. Stability and basin size for the star topology (see fig. [2]). Intermediate delays 
result in large basin volume if averaging is switched on. Shown are the real parts of the 
eigenvalues for the 4 node star motif (dark blue) as well as the basin volume of the producer (dark 
red) and of one consumer (orange) as functions of the delay r for different averaging times: T a = Os, 
Tb = Is, T c = 2s. Parameters a = 0.1/s, K = 8 /s 2 and 7 = 0.25/s were applied. For delay r = 2.1s 
simulations were repeated 21 times, averaged and the standard deviation is shown as a typical error 
bar. 
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with greater delay VbasintV > 0) > Wasin (t = 0) until, for delays r ~ 4s, we obtain close to 
perfect stability with 14asin ~ 1 both for an averaging T b = Is and T c = 2s. In the previous 
paragraph we demonstrated that high averaging times and large delays always destabilize 
the power grid. Hence, we observe a trade-off in curing resonances with averaging and avoid¬ 
ing the rebound effect for delays larger than a critical value r c . Furthermore, basin volume 
reveals that disturbances in a consumer node are less likely to destabilize the system than 
perturbations of the producer (compare dark red and light orange curves in fig. [3]). This 
is intuitively clear as there is only one producer and the topology increases its importance 
even more. 

We conclude that Decentral Smart Grid Control can be applied to the star motif if an 
averaging time of at least T = 2s is used or the price elasticity is smaller than the intrinsic 
damping 7 < a. Additionally, intermediate delays r ~ 4s incorporate the trade-off between 
curing either destabilizing resonances or rebound effects. They increase the basin volume of 
the system and thereby benefit the overall stability of the power grid. 


B. Effect of decentralized production 

In this section we demonstrate that switching from central to decentralized production 
improves the linear stability in the power grid topologies we investigate for small and in¬ 
termediate delays. Specifically, we analyze linear stability for moderately sized lattice and 
cycle networks for different central and decentralized power production. 

For a cycle network decentralization enhances stability significantly (fig. [4]) . For a power 
line coupling of K = 8 /s 2 centralized and decentralized production result in similar stability 
(fig- §3 and e). However, when choosing the critical coupling of the cycle network, i.e., 
the minimal coupling needed so that there exists a fixed point 2 ^, K = 4/s 2 , the cycle with 
central production cannot be stabilized for all considered delays, while this is possible for 
decentralized production (fig. [4]c and f). 

A lattice-like topology for power grids allows stable operation with central power produc¬ 
tion (fig. [ 5 ]). Choosing large couplings of K = 8 /s 2 (fig. |5 ]d and e) or even K = 4/s 2 (fig. 
and g), decentralized and centralized production result in very similar stability. Even when 
operating at the critical coupling of the lattice-like topology K = 2/s 2 , the central power 
production can be stabilized for sufficiently large averaging time T = 4s (fig. [5ji and h). 
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FIG. 4. Central power production in a circle network requires larger capacity K than 
in decentralized power grids. Shown are the ranges of delay r for which the power grid motifs 
with central production (a) or decentralized production (d) are linearly stable. Panels (b) and 
(e) present ranges for a high capacity K = 8 /s 2 , whereas (c) and (f) for K = 4/s 2 . Overall, the 
regions of stability tend to become larger, the larger the average time T. Parameters a = 0.1/s 
and 7 = 0.25/s were applied. 

Note that we chose 7 \ = 0.25/s for all nodes in the networks. Hence, the large producer 
with Piarge = 8 /s 2 adapts relatively less compared to the smaller producers with P sma n = 2/s 2 . 
Nevertheless, the overall adaptation of the whole network is 

N 

AP = ^7HAo;|, (23) 

2=1 

with |Acu| being the maximal angular frequency deviation. Hence, the maximal adaptation 
AP is independent of the power distribution. 

We conclude that a centralized power production requires larger transmission capacities 
compared to a decentralized power production to guarantee stable power grid operation. An 
averaging time of T ~ 4s stabilizes the power grid with Decentral Smart Grid Control for 
all considered delays. Note that our decentralized production utilized short distances to the 
consumers. Decentralized power production tends to allow smaller averaging times, thereby 
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FIG. 5. Central and decentralized power production in a lattice-like topology lead 
to similar stability. In contrast to the cycle network, the lattice-like topology is stable for 
lower coupling K. Shown are the ranges of delay r for which the power grid motifs with central 
production (a) or decentralized production (e) are linearly stable. Panels (b) and (f) present ranges 
for a high capacity I\ = 8 /s 2 , (c) and (g) for an intermediate capacity K = 4/s 2 , finally (d) and 
(h) for K = 2/s 2 . Overall, the regions of stability tend to become larger, the larger the average 
time T. Parameters a = 0.1/s and 7 = 0.25/s were applied. 

offering a greater safe operating space. In addition, a highly connected topology like a lattice 
outperforms the less connected cycle in terms of stability. 

V. SUMMARY AND DISCUSSION 

In this article we applied the concept of ’’Decentral Smart Grid Control” (DSGC), as 
proposed in 2 ^, to different motifs and small networks. We first determined both the linear 
stability and the basin volume of a 4-node-star motif in dependence on the delay time r 
(see equation [ 5 ]) and for fixed averaging times T (see equation [6]). Linear stability analysis 
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reveals two destabilizing effects for the power grid: First, resonance catastrophes destabilize 
the system periodically. This instability can be cured by applying sufficient averaging (fig. 
[2]). Secondly, a rebound effect emerges for large delays and destabilizes the system regardless 
of averaging. The rebound effect sets an upper limit for the delay r = r c and magnitude 
of adaptation response 7 as it has to be smaller than the intrinsic damping of the system 
a. Basin volume analysis gives further probabilistic insight on how well DSGC tames grid 
instabilities. For large averaging times T and intermediate delays r, basin volume approaches 
unity (fig. §■ Hence, for DSGC exists a trade-off in curing resonances with averaging and 
larger delays and avoiding the rebound effect for delays larger than a critical value r c . 

Summarizing the results from linear and basin volume analysis, adaptation has to be 
smaller than the intrinsic damping of the system (7 < a) or the demand response time 
needs to be located in a delay window of safe operation (r < r c ). For values above the 
critical delay r c the system becomes always destabilized, regardless of the averaging time. 
At the same time, averaging and increasing delay is beneficial for system stability in terms 
of basin volume. These results have strong implications on how parameters has to be set for 
real world applications of DSGC. 

In the last section of this article, we demonstrated the usefulness of DSGC with central¬ 
ized as well as with decentralized power production: While it works in both cases, central 
production requires larger line capacities K. For the lattice-like topology, this effect can be 
compensated by using longer averaging times. But decentralized power production is clearly 
advantageous. 

Next research steps include considering heterogeneous networks, i.e. the use of different 
r, 7 , T values for individual nodes, modifying the averaging method, e.g., to a discrete time 
window and extending the DSGC framework to larger network topologies. I 11 this context, 
there remain a couple of open questions that will have to be investigated in more detail, 
namely: What is the reason that we observe delays r for which Fbasin('r) > Wasin^ = 0), 
in particular for larger averaging times? Do we need even larger averaging times T when 
we go to larger networks? How large is the safe operating space to cure instabilities by 
resonances while avoiding the rebound effect for different networks? These are all widely 
open questions. 

In this article, we have demonstrated that Decentral Smart Grid Control constitutes a 
promising control concept, in particular for future power grids that will be more decentralized 
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than the present one. 
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